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GEOMETRIC CONSTRAINED VARIATIONAL CALCULUS. 
I. - PIECEWISE SMOOTH EXTREMALS. 

E. MASSA, D. BRUNO, G.LURIA, AND E. PAGANI 


Abstract. A geometric setup for constrained variational calculus is presented. The analysis deals with 
the study of the extremals of an action functional defined on piecewise differentiable curves, subject to 
differentiable, non-holonomic constraints. Special attention is paid to the tensorial aspects of the theory. 
As far as the kinematical foundations are concerned, a fully covariant scheme is developed through the 
introduction of the concept of infinitesimal control. The standard classification of the extremals into 
normal and abnormal ones is discussed, pointing out the existence of an algebraic algorithm assigning 
to each admissible curve a corresponding abnormality index, related to the co-rank of a suitable linear 
map. Attention is then shifted to the study of the first variation of the action functional. The analysis 
includes a revisitation of Pontryagin’s equations and of the Lagrange multipliers method, as well as a 
reformulation of Pontryagin’s algorithm in hamiltonian terms. The analysis is completed by a general 
result, concerning the existence of finite deformations with fixed endpoints. 


1. Introduction 

Calculus of variations has a very old origin, dating back to the pioneering contributions of Euler, 
Lagrange and Weierstrass * 1 . 

The 20th century witnessed the development of the classical theory towards two different main direc¬ 
tions: on the one hand, the curves on which a given functional is defined, and is required to be minimized, 
were imposed to satisfy a set of differential (or non-holonomic) constraints. This topic took the name 
of constrained calculus of variations. On the other hand, mathematicians accounted for the possibil¬ 
ity, more directly related to optimization theory , of considering also unilateral constraints, described by 
inequalities. 

Both aspects were subsequently merged into control theory. Here, an important role is played by the 
celebrated Pontryagin Maximum Principle, which provides the necessary conditions that must be satisfied 
by a minimizing solution of the given functional. In this connection, besides the original contributions of 
Pontryagin [11], more recent developments have been given by Bellman [12], Hestenes [13], Cesari [14], 
Alekscseev, Tikhomirov and Fomin [15], Schattler and Ledzewicz [16]. 

Recently, Griffiths [17], Hsu [18], Sussmann and his school [19], Sussmann and Liu [20], Mont¬ 
gomery [21], Agrachev and his school [22], Gracia [23], Munoz-Lecanda and his school [24, 25], Jimenez, 
Yoshimura, Ibort, de Leon and Martin de Diego [26, 27, 28] have significantly contributed to develop 
differential geometric approaches to constrained calculus of variations and to control theory. 

This paper deals with constrained variational calculus in the sense defined above. We present a 
fresh geometric approach to the subject, extending known results and suggesting new ones. Attention 
is focussed on the construction of a fully covariant scheme, suited to the study of arbitrary non-linear 
constraints and embodying the possibility of piecewise differentiable extremals. The geometric setup and 
terminology are borrowed from non-holonomic mechanics [29, 30]. 

Throughout the discussion we shall deal with parameterized curves , namely curves whose parameter¬ 
ization is fixed once for all, up to an additive constant; geometrically, this means dealing with sections 
of a fiber bundle, called the event space. In the resulting context, the constraints are described by an 
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Sagan [10]. 
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embedded submanifold of the first jet-bundle of the event space, intuitively representing the totality of 
admissible velocities. 

The study of the extremals is performed using the standard techniques of variational calculus, namely 
discussing the zeroes of the first variation of the action functional with respect to deformations with fixed 
endpoints. Old though the whole argument may seem, some questions are still unsettled: among others, 
the geometrical nature of the concept of normality of a section; the relation between normality and 
deformability; the creation of fully covariant scheme embodying the kinematical and dynamical aspects 
of the theory. 

The arguments are organized as follows: in Sec. 2 the geometric tools needed in the subsequent 
discussion are presented. After a few preliminary remarks, the infinitesimal deformations of an admissible 
section 7 of the event space are discussed via a revisitation of the well known variational equation. The 
novelty consists in the introduction of a suitable transport law for vertical vector fields along 7 , yielding a 
covariant characterization of the “true” degrees of freedom involved in the description of the most general 
admissible infinitesimal deformation. 

The analysis is then extended to arbitrary piecewise differentiable evolutions and to the associated class 
of infinitesimal deformations. The admissible evolutions are classified into ordinary , if every admissible 
infinitesimal deformation vanishing at the endpoints is tangent to some finite deformations with fixed 
endpoints, and exceptional in the opposite case. 

Along the same guidelines, every piecewise differentiable evolution is assigned a corresponding abnor¬ 
mality index , extending and expressing in geometrical terms the traditional attributes of normality and 
abnormality commonly found in the literature for regular evolutions. 

A noteworthy result, established in Appendix A, is the existence of a strict relationship between ab¬ 
normality index and ordinariness, leading to the conclusion that all normal evolutions are automatically 
ordinary. This extends to arbitrary non linear constraints and to piecewise differentiable sections a result 
proved by Hsu [18] in a linear context. 

After these preliminaries, in Sec. 3, attention is focussed on the determination of the broken extremals 
- or extremaloids — of the action functional. Once again, a fully covariant approach is proposed, yielding 
back Pontryagin’s equations [9, 11] as well as the Erdmann-Weierstrass corner conditions [9, 10]. 

The resulting equations are shown to be sufficient for any evolution, and necessary and sufficient for 
an ordinary evolution to be an extremaloid. The same setup is seen to provide a concise approach to the 
Lagrange multipliers method. 

A deeper insight into the nature of the problem is then gained lifting the algorithm to a more natural 
environment, here called the contact bundle , gluing together the velocity space and the phase space. 
This allows recovering the ordinary extremaloids of the original variational problem as projections of the 
extremaloids of a free variational problem on the contact bundle. 

The same environment is also seen to provide an intrinsic characterization of the abnormal evolutions 
in terms of the geometric properties of the contact bundle, as well as a proof of the fact that, under 
suitable regularity assumptions, the original constrained variational problem is locally equivalent to a 
free hamiltonian problem in phase space. 


2. Geometric setup 

In this Section we present a brief review of the geometric tools involved in the subsequent discussion. 
Throughout the paper, we shall freely use the language and methods of differential geometry [31, 32]. 
The terminology will be partly borrowed from classical non holonomic Mechanics [29, 33, 34, 35]. 

2.1. Preliminaries. Let V„+i H denote a fibre bundle over the real line, henceforth called the 
event space , and referred to local fibred coordinates t, q 1 ,..., q n . 

Every section 7 : R — > V n +i is interpreted as the evolution, parameterized in terms of the independent 
variable t, of an abstract system B with a finite number of degrees of freedom. The first jet-bundle 
ji{V n+1 ) —V n +i, referred to local jet-coordinates t,q l ,q l , is called the velocity space. The first jet- 
extension of 7 is denoted by ji( 7 ): M —» ji(V„+i) ■ 
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The presence of differentiable constraints is accounted for by a commutative diagram of the form 

A ji(V„+i) 

( 2 . 1 ) 


Vn+l 


v n+1 


where the submanifold A 


j i(V ri +i) , fibred over V n +i , represents the totality of admissible velocities. 

1 ., z r with transformation laws 

„1 „r\ 


The manifold A is referred to local fibred coordinates t, q 1 ,..., q“, z 
t = t + C, <f =q t (t,q 1 ,...,q n ), Z A = z A {t,q 1 


while the imbedding i : A —> j i(V n +i) is locally expressed as 

(2.2) q l =ip t (t,q 1 ,...,q n ,z 1 ,...,z r ) i = l,...,n, 


with rank 


a ^ 1 •••*»■) 


r. To make it easier, no distinction is made between A and its image 


i{A) C ji(V„+i). 

An evolution 7 : R — > V n +i is called admissible if and only if its first jet extension ^ 1 ( 7 ) factors 
through A, i.e. if and only if there exists a section 7 : R —>• A satisfying j 1 ( 7 ) = i ■ 7 . In the stated 
circumstance, the section 7 , locally described as q 1 = q l (t),z A = z A (t ), is itself called an admissible 
section of A , the admissibility requirement taking the explicit form 


(2.3) 


dq i 

dt 


i> 1 (t,q 1 (t),...,q n (t),z 1 {t),...,z r (t)) . 


With this terminology, the admissible sections of A are in bijective correspondence with the admissible 
evolutions of the system. We shall refer to 7 as to the lift of 7 . 

On account of Eq. (2.3), the admissible evolutions of the system are determined, up to initial data, by 
the knowledge of the functions z A (t). In general, however, these functions have no invariant geometrical 
meaning 2 . 

Perfectly significant tools for the description of possible evolutions of the system are instead the sections 
<7 : V„+i A . Each of them is called a control for the system; the composite map i-a : V n +i —> Ji(V n +i) 
is called an admissible velocity field. In local coordinates we have the representations 

a : z A = z A (t,q 1 ,...,q n ), 

i-a: q l = ip'(t, q 1 ,..., q n , Z A (t, q 1 , , 


pointing out that the knowledge of er does indeed determine the evolution 7 (t) from a given initial 
configuration 7 (^ 0 ) through a well -posed Cauchy problem. 

A control a : V„+i —> A is said to contain a section 7 : R —> V n +i if and only if the lift 7 : H —> A 
factors into 7 = a ■ 7 , i.e. if and only if the jet extension j 1 ( 7 ) coincides with the composite map 
* • er • 7 : R ->• ji(V„+i). 


2.2. Geometry of the velocity space, (i) Given the event space V n +i, we denote by V(V v ^\) the 
vertical bundle associated with the fibration V ra +i —> II and by V*{V n+ 1 ) the corresponding dual bundle. 
The elements of V*{V n + 1 ) are called the virtual 1-forms over V„+i. 

By definition, P*(V„+ 1 ) is canonically isomorphic to the quotient of the cotangent bundle T*(V n +i) 
by the equivalence relation 


(2.4) 


a 


7 r(cr) = n(a f ) 
(j (j oc dt 


The quotient map T*(V n +i) —>• V^*(V n +i) is denoted by w. For each function g G ^(V n +i), the 
section 5g : V n +i —> T^*(V n +i) given by 8g\ x := m(dg\ x ) is called the virtual differential of g. 

Every local coordinate system t,q l in V n +i induces fibred coordinates t,q l ,pi in V^Vn+i), uniquely 
defined by the requirement A = pi(X)Sq^ x VA G Vf (V n +i)• 


2 An obvious exception occurs whenever the manifold A splits into a fibred product V n +i Xr 2, the second factor 
being some fiber bundle over R. 
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Every element belonging to the tensor algebra generated by V(V n +i) and V*(V n + 1 ) is called a virtual 
tensor over V ra+ i. We preserve the notation ( , ) for the pairing between V(V n +i) and V*(V n +i)- 

(ii) The fibred product C(A) := A Xy„ +1 V*(V„+ 1 ) is called the contact bundle: it is at the same 
time a vector bundle over A, isomorphic to the subbundle of T*(A) locally spanned by the 1-forms 
w' := dq l — ijj l dt , and a fibre bundle over V*(V n +i) • The situation is summarized into the commutative 
diagram 

C(A) -—=-»■ V*(V n+1 ) 

(2.5) C | J* 

A —=—>• V n+ i 

We refer C(A) to coordinates t, q\z A ,pi gluing in an obvious way the coordinates in the factor spaces. 
Every a £ C(A) is called a contact 1-form over A. 

As a subbundle of the cotangent space T*(A) , the manifold C(A) is naturally endowed with a Liouville 
1 -form 0 , uniquely defined by the relation 

(Z,0) = ((C k )*(Z),a) V Z £ T a (C(A)) 

and expressed in coordinates as 3 

( 2 . 6 ) 0 = piLO 1 = Pi(dq l — tp l dt). 


(iii) A vector field T over A is called a semispray if and only if its integral curves are jet-extensions of 
admissible sections 7 : R —> V n +i. In coordinates, this property is summarized into the representation 

r = -§- t +r^ + r A 


d 

~dP' 


According to the latter, given any covariant tensor field w on V n +i, the Lie derivative £r( 7 r*(w>)) is 
a covariant field on A, depending only on w and not on the choice of the semispray T. 

The correspondence w (->• £r( 7 r*(w)), essentially identical to Tulczyjew’s operator dx [36, 37, 38], is 
here denoted by and called the symbolic time derivative: by construction, it is a derivation of the 
covariant tensor algebra over V n +i with values in the covariant algebra over A, uniquely characterized 
by the properties 


(2.7a) §-*(.-</)) + -A 

(2.7b) f t (df) = £ r (ir*(<i/)) = dC r (t’(f)) = df. 


(iv) The vertical bundle associated with the fibration A V n +i is denoted by V(A) . Each Y = 
V A (gpr)_ £ 14(A) induces an action on the ring of differentiable functions over V n +i, based on the 
prescription 

e(¥)uy -= Y{f) = Y i^)M)^ v/e ^ (v “ +i) ' 

According to the latter, g(Y) is a derivation of ^(V n +i), he. a (vertical) vector on V n +\- We have 
therefore a fibred homomorphism 

V(A) — 3 —> V (V n+ i) 

1 I 

A —> V n+ l 


expressed in components as 

( 2 . 8 ) 



/ dip 1 \ f & \ 


^For simplicity, we do not distinguish between covariant objects in A and their pull-back in C(A)] namely, we write 
^ for ^ for C*(^) etc. 
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2 . 3 . Vector bundles along sections, (i) Given an admissible section 7: R —> V„+i, we denote by 
V (7) —^ It the bundle of vertical vectors along 7, by V (7) R the restriction of V (.4) to the lift 
7: R — > A , and by A (7) R the bundle of isochronous vectors along 7, meant as vectors annihilating 
the 1-form dt. 

We adopt fibred coordinates t,u l in W(7), t,v A in 1^(7) and t,u\v A in A{ 7), according to the 
prescriptions 


( 2 . 9 a) x = u\X)(-^-\ Vl 6 b( 7 ), 

V^V 7(t( x)) 

( 2 . 9 b) y = v a (Y)(JL'\ Vbebffl, 

V 0 - /y (t(Y)) 

( 2 . 9 c) x = u\X)(-?-\ +V A { X)(°) Viel( 7). 

The symbolic time derivative defined in Sec. 2.2 induces a derivation — still denoted by 4 - of the 

algebra of covariant tensor fields along 7 into the algebra of covariant tensor fields along 7. 

On account of Eqs. ( 2 . 7 a, b), the effect of 4 . on a function f{t) coincides with the ordinary derivative, 
while on a 1 form w = wo(t)(dt)^ + Wi{t)(dq l ) 1 it reads 


( 2 . 10 ) 


dw 

dt 


dwo 

dt 


dt 


dwi 

dt 


dq 1 


d i/d 
~dt 


dt 


dip 1 k dip 1 A 

XT d z 

dq K dz A 


In a similar way, the correspondence ( 2 . 8 ) 
expressed in components as 


determines an injective homomorphism V (7) V (7), 


( 2 . 11 ) 


Q 


dz A ). 


f dip 1 \ / d \ 


(ii) A straightforward argument, left to the reader, ensures that the first jet-bundle j\ (V (7)) is canon¬ 
ically diffeomorphic to the space of isochronous vectors along j 1(7). In jet coordinates, this results into 
the identification 


Z = u\Z) 


_d_ 

~dq l 


u\Z) 


hh)(t(z)) 


d_ 

dtf 


3l (j)(t(Z)) 


vzeMV( 7 )). 


Due to the latter, the push-forward of the imbedding ( 2 . 2 ), restricted to the submanifold A( 7) C 
T(A) , determines a vector bundle homomorphism 


^( 7 ) - - > ji(V( 7 )) 

( 2 . 12 a) 7T* 7T* 

V (7) = V (7) 


expressed in coordinates as 


(2.12b) 



\9z A )~ 


v 


A 


The kernel of the projection A{ 7) —X- V( 7 ) is easily recognized to coincide with the vertical bundle 

V(7). 

(iii) The restriction of the space V*(V n+ 1) to the curve 7 determines a vector bundle V*( 7 ) —A R, 
dual to the vertical bundle V (7). 

The elements of V*( 7 ) are called the virtual 1 -forms along 7. We recall that, by definition, they are 
not 1 -forms in the ordinary sense, but equivalence classes of 1 -forms under the relation ( 2 . 4 ). 

The elements of the tensor algebra generated by 1^(7) and V*( 7 ) are called the virtual tensors 
along 7. 
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Preserving the notation vo : T*(V„+i) —> for the quotient map sending each (ordinary) 1-form 

A along 7 into the corresponding equivalence class and 5q l for the image vo (dq l ), every virtual tensor 
field w : 1R — > V ( 7 ) <8 >]r V* ( 7 ) (g)®, • ■ • is locally represented as 

w = w^.-Xt) 

Every local coordinate system t,q l in V n +i induces fibred coordinates t,pi in V*(j), with 

ft(A):=<(^) 7 , A) VkF*( 7 ). 


2.4. Deformations. A deformation 7 ^ : R —>• V n +i of an admissible section 7 — and, likewise, a defor¬ 
mation 7 £ : 1 R. — ¥ A of the lift 7 — are called admissible if and only if all sections 7 j, 7 ^ are admissible, 
i.e. if and only if they satisfy the condition j 1 ( 7 ^) = * • 75 . 

In coordinates, assuming the representation 7 : q l = q l (t ), = z A (t ), the admissible deformations 

of 7 are described by equations of the form 

7C : = = 

subject to the conditions 

(2.13a) ip l (0 ,f) = q\f), { A (0 ,t) = z A (t), 

(2.13b) ■ 


Setting X l (t) := (^|-) , X A {t) := (^4-) > the infinitesimal deformation tangent to 77 is the 


9 £ ^ ■- V 94 l£ =0 

section X : R —► . 4 ( 7 ) locally expressed as 

1 - xHt) {*\ + xA{t) . 

while the admissibility condition (2.13b) is reflected into the variational equation 


dz A ). 


dX i 


dt 


i _ fdr\ yk , (or\ 

" Wa \dz A )« 


X A . 


(2.14) 

hi. \ t tn •“ ... ........ 

7 v / 7 

The infinitesimal deformation tangent to the projection 7 ^ = 7 r • 7 ^ is similarly identified with the 
section X: R —► P(y) locally expressed by 

d 


X = 7T*X = X l (t) 


dq i 


Collecting all results and recalling Eq. (2.12b) we conclude 


Proposition 2.1. Let 7 : R — > A be the lift of an admissible evolution 7 : R — > V n +i . Then, a section 
X : R — > P( 7 ) represents an admissible infinitesimal deformation of 7 if and only if its first jet extension 
factors through A{ 7 ), i.e. if and only if there exists a section X : R — > A{fj) satisfying ji(X ) = 
Conversely, a section X : R —>• A ( 7 ) represents an admissible infinitesimal deformation of 7 if and only 
if it projects into an admissible infinitesimal deformation of 7 , i.e. if and only if i*X = ji( 7 r*X). 


The proof is entirely straightforward, and is left to the reader. 

Proposition 2.1 points out the completely symmetric roles played by diagram (2.1) in the study of the 
admissible evolutions and by diagram ( 2 . 12 a) in the study of the admissible infinitesimal deformations, 
thus enforcing the viewpoint that the latter context is essentially a “linearized counterpart” of the former 
one. 


Remark 2.1. Throughout the previous discussion, an admissible infinitesimal deformation is meant 
as a vertical vector field along 7 (or as an isochronous vector field along 7 ) tangent to an admissible 
finite deformation whatsoever. A necessary and sufficient condition for this to hold is the validity of the 
variational equation (2.1j). 
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Any additional requirement on the allowed finite deformation — such as e.g. keeping the endpoints 
fixed — is reflected into a corresponding restriction on the infinitesimal ones. We shall return on this 
point in Sec. 2. 7. 


2.5. Infinitesimal controls, (i) According to Proposition 2.1, the admissible infinitesimal deformations 
of an admissible section 7 : R —> V n +i are in 1 1 correspondence with the sections X : R — > A ( 7 ) 
satisfying the consistency requirement i*X = ji(n*X). 

In local coordinates, setting X = X l {t) + X A (t) , the requirement is expressed by the varia¬ 
tional equation (2.14). 

Exactly as it happened with Eq. (2.3), Eq. (2.14) indicates that, for each admissible X : the functions 
X l (t) are determined, up to initial data, by the knowledge of X A (t). Once again, however, one has to 
cope with the fact that the components X A have no invariant geometrical meaning. 

The difficulty is overcome introducing a linearized version of the idea of control. 

Definition 2.1. Let 7 : 1R, —> V n +i denote an admissible evolution. Then: 

• a linear section h : V ( 7 ) —> A ( 7 ), meant as a vector bundle homomorphism satisfying 7r* -h = id, 
is called an infinitesimal control along 7 ; 

• the image TL( 7 ) := h(V("f)), viewed as a vector subbundle of A ( 7 ) —> R, is called the horizontal 
distribution along 7 induced by h; every section X : 1R, —> A ( 7 ) satisfying X(t) £TL ( 7 ) V t £ It 
is called a horizontal section. 

Remark 2.2. The term infinitesimal control is intuitively clear: given an admissible section 7 , let o : 
Vn+i —s- A denote any control containing 7 , i.e. satisfying ( 7-7 = 7 . Then, on account of the identity 7 r* • 
cr* = (77 ■ a) * = id , the tangent map o * : T (V n + 1 ) —>■ T(A) , restricted to V ( 7 ), determines a linear section 
< 7 * : V ( 7 ) —>• A ( 7 ) satisfying the requirements of Definition 2.1. The infinitesimal controls may therefore 
be thought of as equivalence classes of ordinary controls having a first order contact along 7 . 


(ii) Given an infinitesimal control h : V(j) —> A ( 7 ), the horizontal distribution TLifi) and the vertical 
subbundle V ( 7 ) split the vector bundle A ( 7 ) into the fibred direct sum 

(2.15) A{fl) = H(7)® r V(7), 

thereby giving rise to a couple of homomorphisms Vu '■ A ( 7 ) —> TLifl) (horizontal projection) and 
Vv : A( 7 ) —^► V( 7 ) (vertical projection), uniquely defined by the relations 

(2.16) Vu — h ■ 7r* ; Vv = id — Vu ■ 

In fibred coordinates, every infinitesimal control h : V ( 7 ) —> A ( 7 ) is locally represented as 

(2.17) v A = hi A (t)u l . 


In this way: 

• the horizontal distribution LI ( 7 ) is locally spanned by the vector fields 


(2.18) 


di := h 


_d_ 

'dtp 


_d_ 

~dq l 




(2.19) 


7 \ / 7 

• every vertical vector field X = X l (t)(-g~) along 7 may be lifted to a horizontal field h(X) 
along 7 , expressed in components as 

d\ A f d \ 


h{X) = X\t)di = X\t) 


dq l 


\dz A J. 


• every vector X = X z (-£j T )„ + A' A (g|x). £ A(fi) admits a unique representation of the form 
X = Vu{X)+V v (X), with 

( 2 . 20 ) r n (x) = x‘d u v v{ x ) = (x A -xnA)(Xy =:U *(J pi y; 
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• every 1 -form v along 7 determines a linear functional on . 4 .( 7 ) and therefore, by duality, a 
virtual 1 -form h*(v) along 7 , uniquely defined by the relation 


(2.21) (h*{u),X) = (u,h(X)) Vie V( 7 ) =► h*(v) = (v, d { ) Sq ^. 


The role of Definition 2.1 is further enhanced by the following 

Definition 2.2. Let h be an infinitesimal control along the (admissible) section 7 . A section 
X : R — y V( 7 ) is said to be h-transported along 7 if and only if its horizontal lift h(X) : R — > A(ff) is 
an admissible infinitesimal deformation of 7 , i.e. if and only if i * ■ h(X) = ji(X). 

In coordinates, setting X = X l (t)(-^ T ) and recalling Eqs. (2.14), (2.19), the condition for h- 

transport is expressed by the linear system of ordinary differential equations 

X k = X k d k ip i . 

From the latter, on account of Cauchy theorem, we conclude that the /i-transported sections of V ( 7 ) 
form an ?r-dimensional vector space V),, isomorphic to each fibre ( 7 )|* through the evaluation map 
X X(t). We can therefore state: 

Proposition 2.2. Every infinitesimal control h : V(pf) —> A{ 7 ) determines a trivialization of the vector 
bundle P( 7 ) U R. 


( 2 . 22 ) 


dX l 

dt 


dip 1 

dq k 


TTT ) + hk A 


I di/’*! 

\d^L 


Proposition 2.2 provides an identification between sections X : R —» V("f) and vector valued functions 
X : R — > Vh — whence, by duality, also an identification between sections A : R —>• W(7) and vector 
valued functions A : R —> Vj* — thus opening the door to the introduction of an absolute time derivative 
for virtual tensor fields along 7 . The algorithm is readily implemented in components. 

To this end, let {e( a )}, {e^ } be any pair of dual bases for the spaces 14, Vf. By definition, each e( Q ) is 
a vertical vector field along 7 , obeying the transport law (2.22). In coordinates, setting e( a ) = e fa (■£■?) , 
this implies the relation 


(2.23a) 


de,\ 

(a) 

dt 


dkff ■ 


In a similar way, each is a virtual 1-form along 7 , expressed on the basis 5q l i 7 as e ^ 5q l | 7 , 

with e < 'f l ed-j = <5^ . On account of Eq. (2.23a), the components e/”' 1 obey the transport law 


(2.23b) 




= 0 


def 


(a) 


dt 


= -e^f diip 0 . 


The functions 
(2.24a) 


de {a) ■ 
dt (“) _ 



de,\ 

(a) 

dt 


are called the temporal connection coefficients associated with the infinitesimal control h in the coordinate 
system t,q l . Eqs. (2.18), (2.23a,b) provide the expression 


(2.24b) 


rs.-w.-m) 


\ d( i l A 


\dz A J~ ' 


Given any section X : R —>• V ( 7 ), the definition of the operator A yields the evaluation 


DX d/ i a \ 

= — ( X, e ( ’ ) e( Q ) = 


Dt 


dt 




written more simply as 


(2.25a) 


DX 

~Dt 



+ rr, 3 j 
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In a similar way, for any virtual 1 -form A : 1R, —> V* ( 7 ), the same argument provides the result 


(2.25b) 


DA d . , 

— = — ( A, e< 


Dt 


dt 


e (a) ) e(a) = ^( A i e (a)) e ? V | 7 = (^ - V| 7 - 


(iii) In view of Eq. (2.20), every infinitesimal deformation X of the section 7 admits a unique represen¬ 
tation of the form X = h(X) + U, where: 

• X = tv*X := X 1 (^) is a vertical field along 7, namely the (unique) infinitesimal deformation 
of 7 lifting to X ; 

• U = Tv{X) := D A (g^-). is a vertical vector field along 7. 

In terms of this decomposition, the variational equation (2.14) reads 


dX l 

dt 


= X k 


dtp* 

dq k 


(IS) w + v*). 


\dz A )- 


7 v '7 

Recalling Eqs. (2.24b), (2.25a), as well as the representation (2.11) for the homomorphism 
q : V (7) —> V (7), the latter equation is more conveniently written as 

f) V 

(2.26a) 


Dt 


= e{U) = g{Vv(X)) 


or also, setting X = X a e ^ a ), U = U A (^73-) - and expressing everything in the h-transported basis e^ a 

dX ‘ 


(2.26b) 


= (e (a) , 


dt \“ ’ KV “' / / \dZ A J-y 

Precisely as the original equation (2.14), Eq. (2.26a) points out that every infinitesimal deformation 
X is determined by the knowledge of a vertical vector field U : R —» V (7) through the solution of a well 
posed Cauchy problem: the advantage of the present formulation is that all quantities have now a precise 
geometrical meaning relative to the horizontal distribution TL (7) induced by the infinitesimal control h. 

As we shall see, this aspect has useful consequences. At the same time, of course, one cannot 
overlook the fact that, in the formulation of a typical variational problem, no distinguished section 
h : V (7) — > A (7) is generally included among the data, and no one is needed in order to formulate the 
results. 

In this respect, the infinitesimal controls play the role of gauge fields, useful for covariance purposes, 
but unaffecting the evaluation of the extremals. Accordingly, in the subsequent analysis we shall regard h 
as a user-defined object, eventually checking the invariance of the results under arbitrary changes h —>• h'. 


2.6. Corners. As anticipated in the Introduction, in the study of the extremals of a variational problem 
we shall not only consider ordinary sections, but also piecewise differentiable ones, defined on closed 
intervals. To this end, we stick to the following standard terminology: 

• an admissible closed arc (7, [a, b]) in V n +i is the restriction to a closed interval [a, b] of an 
admissible section 7 : (c, d) —> V n +i defined on some open interval (c, d) D [a, b ]; 

• a piecewise differentiable evolution of the system in the interval [to,t\] is a finite collection 

( 7 , [*o,*i]) : = {(l (s \ [«s-i,«s]), S=1,...,N, t 0 = a 0 <a x < ■ ■ ■ < a N = h} 
of admissible closed arcs satisfying the matching conditions 
(2.27) 7 (s) W = 7 (s + 1) M Vs = l,...,iV — 1. 

According to Eq. (2.27), the image 7(t) is well defined and continuous for all to <t <t\. This allows 
to regard the map 7 : [to,D] —> V n +i as a section in a broad sense. The points 7(to), 7(ti) are called 
the endpoints of 7. The points x s := 7 (a s ), s = 1,. .., N — 1 are called the corners. 

Consistently with the stated definitions, the lift of an admissible closed arc (7, [a, b ]) is the restriction 
to [a, 5] of the lift 7 : (c,d) —> A, while the lift (7, [to,ii]) of a piecewise differentiable evolution 
{(7 (s) , [a s _i,o s ]) } is the family of lifts 7^, each restricted to the interval [a s _i,a s ]. 
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With this definition, the image 7 (t) is well defined for all t ^ a 1,... thus allowing to regard 

7 : [to,ti] —> A as a (generally discontinuous) section of the velocity space. In particular, due to the fact 
that the map i : A — > ji(V n + 1) is an imbedding of A into an affine bundle over V n +i, each difference 

[t] 0s = *(7 (s + 1) 0s)) - *(7 (s) (a s )), s = 1 

identifies a vertical vector in T Xs {V n+ 1), called the jump of 7 at the corner x s . 

In local coordinates, setting gU(f) '■= Q z ( 7^(t)) and introducing the notation 
[^(t)] a := V ,I (7^ S + 1 H a s)) ~ V’ l (7^ s H a s)) to indicate the jump of the function i/> z (y(t)) at t = a s , 
Eqs. (2.3), (2.27) provide the representation 


(2.28) 



(m.-mM) 



Along the same guidelines, an admissible deformation of an admissible closed arc (7, [a, 6]) is a 
1-parameter family (7^, [a(£), 6(£)]), |£| < e of admissible closed arcs depending differentiably on £ and 
satisfying (70, [a(0), 6(0)]) = (7, [a, 6]). Notice that the definition explicitly includes possible variations 
of the reference interval [a(£), 6(£)]. 

In a similar way, an admissible deformation of a piecewise differentiable evolution (7, [to, ti ]) is a col¬ 
lection {(7^, [a s _i(£),a s (£)])} of deformations of the various arcs, satisfying the matching conditions 

(2.29) 7 ( / } M0) = 7 ( / + 1) M0) V |£| < e, s = 1,... ,!V — 1. 

In the stated circumstance, the lifts 75 and 7^, respectively restricted to the intervals [a(£),6(£)j 
and [a s _i(£), a s (£)) are easily recognized to provide deformations for the lifts 7 : [a, 6] —» A and 
7 (s) : [o s _i,a s ] -» A. 

Unless otherwise stated, we shall only consider deformations leaving the interval [fo,ti] fixed, i.e. sat¬ 
isfying ao (£) = to , in (£) = *i V £ . No restriction will be posed on the functions a s (£), s = 1,..., N — 1. 
Each curve x s (£) := 7^(0*, (£)) will be called the orbit of the corner x s under the given deformation. 

In local coordinates, setting = pX(£,t ), the matching conditions (2.29) read 

(2-30) ¥>&)(£, «»(£)) = Vl + i)(£, «»(£)), 

while the representation of the orbit x s (£) takes the form 

(2.31) x s (£) : t = a s (£), q l = V> («)(£> ««(£))• 


The previous arguments have a natural “infinitesimal” counterpart. More specifically: an admissible 
infinitesimal deformation of an admissible closed arc (7, [a, 6]) is a triple (a,X,/3), where X is the 
restriction to [a, 6] of an admissible infinitesimal deformation of 7 : (c, d) —> Vn+i, while a, j3 are the 
derivatives 


(2.32) 


a = 


da 

d£ 


P = § 

£=0 ^ 


£=o 


expressing the speed of variation of the interval [a(£), 6(C)] at £ = 0. 

In a similar way, an admissible infinitesimal deformation of a piecewise differentiable evolution (7, [to, f 1]) 
is a collection | • a s _i, A'( s ), a s ■ ■ • } of admissible infinitesimal deformations of each single closed arc, 

with a s = ^-| „_ 0 and, in particular, with a 0 = &n = 0 whenever the interval [to,ti] is held fixed. 

Let us analyse the situation in closer detail. To start with, we notice that the quantities a s , X( s ) are 
not independent: Eqs. (2.30) imply in fact the identities 






d( P l (s) das_ 

dtj 


dip 


is) 


d^ls + i) das_ 


a£ dt < a£ dt 

From these, evaluating everything at £ = 0 and recalling the relation between finite deformations and 
infinitesimal ones, we get the jump relations 


(2.33) 


X 


(® + i) 


- X 


w)c 


= -a s 


dq 


(« + l) 


dq 


is) 


dt 


dt 


= ~a s 


V’Xt) 
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Furthermore, the admissibility of each single infinitesimal deformation X( s ) requires the existence of 
a corresponding lift X(s) = Xfa (^r)^(s) + (jpx) - ( s ) satisfying the variational equation (2.14). 

Both aspects are conveniently accounted for by assigning to each 7 an (arbitrarily chosen) infinites¬ 
imal control h^ : ^(7^) —> 4(77. In this way, proceeding as in Sec. 2.5 and denoting by (-^) 7 («) 
the absolute time derivative along 7 ^ induced by h^ we get the following 


Proposition 2.3. Every admissible infinitesimal deformation of an admissible evolution ( 7 , [to,ti]) 
over a fixed interval [io,ii] is determined, up to initial data, by a collection of vertical vector fields 
{ t/( s ) = Ufy (gfjr)^(s) }j s = 1,..., N and by N — 1 real numbers oq, .. . , apf-i through the covariant 
variational equations 


(2.34) 


DX 


M 


Dt 




d_ 

<y(») W 


s = l,...,N 


■ .y \“/ \ / -y V 3 / \ " '1 / -y (®) 

completed with the jump conditions (2.33). The lift of the deformation is described by the family of vector 
fields 


(2.35) 


X {s) = h^(X {s) )+ U {S) , s = 


The proof is entirely straightforward, and is left to the reader. Introducing n piecewise differentiable 
vector fields di,...,d n along 7 according to the prescription 

dfit) = h {s) (J^r) Vte(a s _i,a s ), s = l,...,N, 

Eq. (2.35) takes the explicit form 

(2.36) X M = h<-> (A) J + U M = JC ( i, a. + Ufa (^)_ 

on each open arc 7^' : (a s -i,a s ) —>• A. 

To discuss the implications of Eq. (2.34) we resume the notation ^(7) and respectively for 

the totality of vertical vectors and for the totality of virtual 1-forms along 7 4 . We then define a 
transport law for virtual tensors on 7, henceforth called /(-transport, gluing h ( s )-transport along each 
arc (7^, [a s _i,a s ]) and continuity at the corners, i.e. continuity of the components at t = a s . 

In view of Proposition 2.2, the /(-transported vector fields form an n-dimensional vector space 14 , 
isomorphic to each fibre W (7) |* • Likewise, the /(.-transported virtual 1-forms span a vector space V)*, 
dual to 14 and isomorphic to V/. 

This provides canonical identifications of V (7) with the cartesian product [to, t \] X V), and of V* (7) 
with the product [to,ti] x V)*, thus allowing to regard every section X : [to,ti] —> V{j) as a vector 
valued map X : [t 0 ,ti] —1 14 , and every section A : [t 0 ,ti] —> as a map A : [to,ti] —> Vjf. 

Exactly as in Sec. 2.3, the situation is formalized referring 14 and Vf to dual bases {e( a )}, {e^} 
related to (g|r) > Sq l | 7 by the transformation 

(2.37) = e ( “ 3 (t) e (o) , <V| 7 = e ( * } (t) e (a) . 

Given any admissible infinitesimal deformation {(X( s ),ck s )}, we now glue all sections 
X( s ) : [a s _i,a s ] —> V(7^) into a single, piecewise differentiable map X : [to,ti] —> 14, with jump 
discontinuities at t = a s satisfying Eq. (2.33). 

For each s = 1 ,... ,N we have then the relations 

(2-38) X (s) = X a (t)e {a) , (s) = Vte (a s _i,a s ). 

In a similar way, we collect all fields £/( s ) into a single object J7, conventionally called a vertical vector 
field along 7. 


^Notice that this definition is perfectly meaningful even at the corners 7 (a s ). 
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In this way, the covariant variational equation (2.34) takes the form 


(2.39a) 


dX 


dt 


- = cV“4f4) 


\dz A )~ 


completed with the jump conditions 

(2.39b) [X a ] a = ef(a a ) = -a s ef{a s ) [^(7) 


s = 1. 


2.7. The variational setup, (i) Given an admissible piecewise differentiable section 
7 : [to, ii] —> V„ + i, let SO denote the infinite-dimensional vector space formed by the totality of vertical 
vector fields U = , s = 1,..., AT} along 7. Also, let 2H denote the direct sum SO © R N_1 . 

By Eqs. (2.39a,b), every admissible infinitesimal deformation A of 7 is then determined, up to initial 
data, by an element ( U , a) := (£/, aq,..., ajv-i) G 2U. 

For variational purposes, let us now focus on the class of infinitesimal deformations X : [to, ti] —> V (7) 
vanishing at the endpoints of'). Setting X(to) = 0, Eqs. (2.39a,b) provide the evaluation 


(2.40) 


X(t) = 




E 

a s <t 


( a ) / 

e- i { a s 


^(7) 


^(a 


Denoting by T : 22J —> Vh linear map defined by the equation 


(2.41) 


T {U,a) = {^J t uA ^i ] [^Zx) dt - as e( * 0) ( as ) ^( 7 ) 


N-l 


e (a) j 


it is then an easy matter to conclude that the required class is in 1 1 correspondence with the subspace 
ker(T) C 2U. 

Actually, as anticipated in Remark 2.1, what really matters in a variational context is not the space 
ker(T) itself, but the (possibly smaller) subset X C ker(T) formed by the infinitesimal deformations 
tangent to admissible finite deformations with fixed endpoints. 

The linear span of X, henceforth denoted by A (7), will be called the variational space of 7. The 
evolutions of the system will be classified into ordinary , when A(y) = ker(T) and exceptional , when 
A (7) C ker(T). 


(ii) A further important classification comes from the nature of the inclusion T (2U) C 14: an evolution 
(7, [to,ti]) is called normal if T(2U) = 14, abnormal in the opposite case 0 . It is called locally normal 
if its restriction to any closed subinterval [t' Q . t ,\] C [f 0 ,4] is normal. The dimension of the annihilator 
(Y(2U))° C 14* is called the abnormality index of 7. 

In this connection, a useful result is provided by the following 


Proposition 2.4. The annihilator (Y(2U))° C 14* coincides with the totality of h~transported virtual 
1-forms p = pi Sq l \-f satisfying the conditions 


(2.42a) = 

(2.42b) p l {a s )[ip l { 7)] aa = 0, 


A = 1,... ,r, 
s = 1,... ,iV — 1. 


Proof. In view of Eq. (2.41), the subspace (Y(2U))° C 14* consists of the totality of elements p = 
p a e = p a e^ Sq l | T satisfying the relation 





N-l 


dt - 


e ( f } (« s ) W{t) 


= 0 


v(£/,a)e an, 


clearly equivalent to Eqs. (2.42a,b). 


□ 


^ As we shall see, when applied to the extremals of an action functional, the terminology agrees with the current one 
(see, among others, [21, 20] and references therein). 
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On account of Eqs. (2.24b), (2.25b), the condition of /i-transport of p along each arc 7W is expressed 
in coordinates as 

' dip k \ , u a f 

dz A 


(2.43) 


dpi , 

1T + Pk 


dq i 


= 0, 


hi Pk 

■ 7 \ / 7 

the cancellation being due to Eq. (2.42a). The content of Proposition 2.4 is therefore independent of the 
choice of the infinitesimal controls h ^ : V( 7^) —> A (7^'). 

This aspect is formalized denoting by 7 = (gj) + ip '.* (g|») 7 the tangent vector to the curve 7. 
Recalling the expression (2.10) for the symbolic time derivative of a 1-form along 7, we have then the 
following 

Corollary 2.1. The annihilator (Y(2R)) 0 is isomorphic to the vector space formed by the totality of 
continuous (ordinary) 1-forms p = po(t)dt^ + Pi(t)dql along 7 satisfying the conditions 

dp 


<'•+>-»• !=»■ 


(2.44) 

Proof. In coordinates, Eqs. (2.44) are summarized into the system 
po = ~ Pi ip |\ , 


0 = ~T f (pPP&dttf + ^ dqfc + Pl dipfo = 


dt 

dpi 

dt 


Pk 


dip 1 ' 

~dq i 


It 


(dq l - ip'dtj + p k 


dip k \ 


A dz A 

dz A -— dt 

dt 


If 


From the latter it is easily seen that the quotient map vj : T*(V n +i) V*{^) sets up a 1-1 correspon¬ 
dence between continuous 1 form satisfying Eqs. (2.44) and virtual 1-forms p := pi(t)Sq 1 ^ satisfying 
Eqs. (2.42a,b), (2.43). ’ □ 

Straightforward consequences of Proposition 2.4 are: 

• in the absence of constraints, all evolutions are automatically normal ; 

• the abnormality index of a piecewise differentiable section 7 cannot exceed the abnormality index 
of each single arc 7^ . 

These and other topics related to normality are discussed in Appendix B. 

An important insight into the hierarchy between normality and ordinariness is provided by the fol¬ 
lowing 

Proposition 2.5. The normal evolutions form a subset of the ordinary ones. 

The result is proved in Appendix A. In this connection, see also [18]. 


3. Calculus of variations 

3.1. Extremaloids. Let us now come to the central problem of variational calculus. Let «S? € ^(A) 
denote a differentiable function on the velocity space A, henceforth called the Lagrangian. Also, let 
(7, [to, *i]) (7 for short) denote an admissible piecewise differentiable evolution of the system, defined 
on a closed interval [to,ti] C R. 

Indicating by 7 the lift of 7, consider the action functional 

(3.1) J[ 7 ] := [s?dt := J2 r <£)dt 

J 7 s—l Ja.a-1 

Definition 3.1. The evolution 7 is called an extremaloid for the functional (3.1) if and only if, for 
every admissible deformation 7^ = {(7^, [a s _i(£), a s (£)]) } with fixed endpoints, the function 

1(0 := / J?dt = J2 / (7 fY(X)dt 

Jit s= i i(f) 


is stationary at ( = 0. 
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The content of Definition 3.1 is formalized denoting by X( s ) the infinitesimal deformation associated 


~ (s) rl 

with each 7^ • Recalling Eq. (2.36) as well as the definition a s = , we have then the evaluation 


da s 


(3.2a) 


dl 


N 

= E 

S=0 s=l 


a(£) 


J a s _i(£) 




€=o 


N 


= E 


X( s) {^)dt + a s ^f(7 (s) (a s )) - a s _i^f(7 (s) (a s _i)) 

s=l v J “*-1 

On account of the requirement ao = &n = 0, denoting by 

[^(7)] aa := [J5f(7 (s + 1} (« s ))-^(7 (s) (« s ))] 
the jump of the function _S?(7(f)) at t = a s , Eq. (3.2a) may be concisely written as 
dT N /■“» / f) \ iv ~ 1 

< 3 - 2b > W =E/ 4)W + ! '(.,a ? 

J=0 8=1 “ a s — l OZ g=1 

To handle Eq. (3.2b), we introduce AT virtual 1 forms = p^(t)Sq l \^ (s) (one for each arc 7^) 
satisfying the transport law 6 

Dp is) ^ = h*(dJ?) Ms) = (a*JSf). (s) <V| 7 w 


(3.3a) . „ , 

v Dt ; 7 (8) 

as well as the matching conditions 


(3.3b) 


j(«) I = pO + i) 


s = 1,..., TV — 1. 


Once again, for notational convenience, we collect all p ^ into a continuous, piecewise differentiable 
section p : [t 0 ,ti] —> according to the prescription 

(3.4) p(t) = p^\t) V t G [og_i, a g ]. 

On account of Eqs. (3.3 a,b), p is then uniquely determined by 2z?, up to an arbitrary Ti-transported 
virtual 1-form along 7. 

Taking the covariant variational equation (2.34) as well as the duality relations 
(^t) ^ , < 5<7 fc | 7 o)^ = S y - into account, by Eq. (3.3a) we get the expression 


xl s) d i X=(X 


■is), 


Dp^ 

Dt Us)/ dt 


X( s ), P (s) ) - 


DX 


(s) 


Dt Us) 


i(») \ - 


whence also 


dt 


x/ s) 8i(&) dt = 


X is) Pi 


is) 


(X i r>^\ — TjA 

\ X is) Pi J p, \dz A ) Al is) {s) 


Summing over s and recalling Eqs. (2.33), (3.3b) as well as the conditions X(t-o) = X(ti) = 0, this 
implies the relation 

N po-s _ rti / r 

E / x is) dt = - p t ( -E ) U A dt+J2a s ^(7) Pi{a s ). 

g=l " &s —1 Jt 0 \OZ J 7 S=1 L Ja * 

In this way, omitting all unnecessary subscripts, Eq. (3.2b) attains the final form 

nt\ / r\ r/> c\ I i \ N — l 


(3.5) 


dl 

d£ 


£=o Jt Q 


J (§pr - P* ^4) uA dt+ E “s [piiV^i 7) - -^(7) 


® For the notation, see Eq. (2.21) 
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In the algebraic environment introduced in Sec. 2.7, the previous discussion is naturally formalized 
regarding the right-hand side of Eq. (3.5) as a linear functional dX 1 : 22J —> R on the vector space 
211 = QJ©R JV - 1 . 

A necessary and sufficient condition for 7 to be an extremaloid for the functional (3.1) is then the 
vanishing of dl 1 on the subset X C 2U formed by the totality of elements (Z7, a) tangent to admissible 
finite deformations with fixed endpoints. By linearity, this condition is mathematically equivalent to the 
requirement 

(3.6) A(y) C ker (dX 7 ), 

A (7) = Span(X) C ker(T) denoting the variational space of 7. 

As we shall see, Eq. (3.6) provides an algorithm for the determination of all the extremaloids of the 
functional (3.1) within the class of ordinary evolutions. 

The exceptional case is more complicated: the lack of an explicit characterization of the space A (7) 
in terms of local properties of the section 7 requires in fact a direct check of Eq. (3.6) on each exceptional 
evolution. 

In what follows we pursue an intermediate strategy; namely, rather than dealing with Eq. (3.6) we 
discuss the implications of the stronger requirement 

(3.7a) ker(T) C ker(dZy). 

According to the classification introduced in Sec. 2.7, the latter is necessary and sufficient for an 
ordinary evolution, but merely sufficient for an exceptional one, to be an extremaloid of the functional 
(3.1). 

By elementary algebra, the requirement (3.7a) is equivalent to the existence of a (possibly non-unique) 
linear functional <7 : —> R — i.e. of a /i-transported virtual 1 form along 7 — satisfying the relation 

T 

2U -- V h 

a 
R 



Setting a = a a e ^ = cq 5q l | 7 and recalling Eqs. (2.41), (3.5), Eq. (3.7b) amounts to the condition 


’to 


f d££ \ rrA 

V dz A Pi dz A ) 


N -1 


S— 1 


U A dt + ^2 a s Pi(t)r( 7) - Af(q) 


U 


'to 




dt - ^2 a s e2\a s ) 

f dr \ 


T S =1 

t\ 


N-l 


'to 


UAai [j)ZAj„ dt - a s^(a s ) n 7 ) 

By the arbitrariness of (£/, a ), the latter splits into the system 
fit/* Bib 1 

( 3 - 8a ) —-(Pt+vb/-. 4=0, A = l,...,r, 


(3.8b) 


dz A v dz A 

(.Pi +&i)r( 7) - Af(y) 


= 0. 


s = l,...,N -1. 


Collecting all results, and recalling Propositions 2.4, 2.5 we conclude 

Theorem 3.1. Given an admissible evolution 7, let p( 7) denote the totality of piecewise differentiable 
virtual 1-forms p = pt(t) Sq l | 7 along 7 satisfying Eqs. (3.3a,b), (3.4) as well as the finite relations 


(3.9a) 


Pi 


dr cCzf 
dz A dz A ’ 


A = 1,.... 
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and the matching conditions 


(3.9b) 

Then: 


Pi^ l { 7) - -^( 7 ) 


= 0, 


s = 1,..., TV — 1. 


a) the condition p(y) 7 ^ 0 is sufficient /or 7 to 6 e an extremaloid for the functional (3.1); 

b) if 7 is an ordinary evolution, the condition p(y) 7 ^ 0 is also necessary for 7 to be an extremaloid; 

c) 7 is a normal extremaloid, namely an extremaloid belonging to the class of normal evolutions, if 
and only if the set p(y) consists of a single element. 


Proof. In view of Eqs. (3.5), (3.9a,b), the existence of at least one p € p(y) ensures the validity of 
^|^_ 0 = 0 for all admissible infinitesimal deformations vanishing at the endpoints of 7 . Assertion a) is 
then a direct consequence of Definition 3.1. 

In particular, according to our previous discussion, if 7 is an ordinary extremaloid, for any continuous 
virtual 1-form p = Pi5q l i 7 obeying the transport law (3.3a) there exists at least one h-transported 
1-form <7 = <Ji 5q z \j satisfying Eqs. (3.8a,b). The sum p + a = (pi + cq) Sq is then automatically in 
the class p(y), thus proving assertion b). 

Finally, as pointed out in Sec. 2.7, the normal evolutions form a subclass of the ordinary ones, uniquely 
characterized by the requirement (Y(2D))° = {0}. 

Therefore, according to assertion b), a normal evolution 7 is an extremaloid if and only if the 
class p(y) is nonempty. Moreover, by Eqs. (3.3a), (3.8a), given any pair p,p' G p(y), the differ¬ 
ence p — p' is an h-transported 1-form satisfying Eqs. (2.42a,b). By Proposition 2.4 this implies 
p — p' € (T (2D )) 0 => p = p ', thus establishing assertion c). □ 


In view of Eqs. (2.24b), (2.25b), for any p € p( 7 ) the transport law (3.3a) simplifies to 

'dif k ' 


dpi 

dt 


Pk 


dq i 


+ hi A pk 



(d&\ 

= \WJ* 



7 ^ ' '7 ' ^ '7 

the cancellation being due to Eq. (3.9a). Exactly as it happened with Proposition 2.4, all assertions of 
Theorem 3.1 have therefore an intrinsic meaning, irrespective of the choice of the infinitesimal controls 

ftW : ^(yW)->• A(7W). 

The previous arguments show that the determination of the ordinary extremaloids of the functional 
(3.1) relies on 2 n + r equations 

dq i 
dt 

dpi 


(3.10a) 

(3.10b) 

(3.10c) 


= il>\t,q\z A ), 


dip k _ dT£ 

dt + ~W Pk ~ 


Pz ■ 


dif* 


dz A dz A 

for the unknowns q l (t) , Pi(t ), z A (t), completed with the continuity requirements 

(3.11) |y] n . = [p ,] a3 = [PiV-X]*, =0, s = 1,... ,N — 1. 

As already pointed out, all equations are independent of the choice of the infinitesimal controls, and 
involve only the “true” data of the problem, namely the Lagrangian Jzf and the constraint equations 
(2.2). In particular: 

• the algorithm (3.10a,b,c), (3.11) is invariant under arbitrary transformations of the form 


(3.12) 


JSf -»■ JSf 


df 


d f 


ib K 

dt dq k ^ ’ 


Piit) -> Pi(t ) 


df\ 

dq 1 ) 


7 (t) 


f(t , q 1 ,..., q n ) being any differentiable function over V n +i; 

the last pair of equations (3.11) extend to the present context the well known Erdmann- Weierstrass 
corner conditions of holonomic variational calculus [10, 9]. 
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In strict analogy with the result established in Corollary 2.1, a more compact description of the 
extremality conditions is provided by the following 

Corollary 3.1. The class p(j) is in 1-1 correspondence with the (possibly empty) affine space p*(y) 
formed by the totality of continuous 1-forms p = p^dt^ + Pidql along 7 satisfying the conditions 

dp 


(3.13) 


(p, 7 ) = ■% > 


dt 


= CLSfj^. 


Proof. In coordinates, Eqs. (3.13) are summarized into the system 
(3.14 a) p 0 = .5% - pi if fa, 

0 = (■% - Pii>fo) <% - ^ dqfa - pi dtpfa = 


dt 


dq i )\i 


(3.14b) 


(d&\ dpi fdif k 

A )i dt Pk \t 

(d&\ _ (dif k \ 
\dz A ) I7 Pk \dz A ) 17 


(dq l — 'if' 1 dt^j 


14 


A dz A 

dz A - dt 

dt 


14 


Eq. (3.14a,b), together with the assumed continuity of the components po ,pi at the corners, reproduce 
the content of Eqs. (3.10b,c), (3.11). Recalling the previous discussion on the characterization of the 
class p(y), we conclude that the quotient map w : T*(V n + 1 ) —> 14* ( 7 ) sets up a 1-1 correspondence 
p’(7)«p(i □ 


In view of Corollary 3.1, all assertions of Theorem 3.1 may be rephrased, systematically replacing the 
class p(q) with p*(y). 

For future reference we point out that, on account of Eqs. (3.10b,c), (3.14a), the component Po(t) 
satisfies the evolution equation 


(3.15) 


dpp _ dJftf _ dp, j _ #17 

dt dt dt ^ ^ 1 dt 


formally similar to Eq. (3.10b). 




3.2. Lagrange multipliers. As an illustration of the algorithm developed in Sec. 3.1, we now discuss a 
constructive approach to the Lagrange multipliers method. To this end, we consider a situation in which: 

• the action integral involves a “free” Lagrangian C, viewed as a function on ji(V n -|_i), with local 
expression C = C(t, q l , q l ); 

• the presence of the constraints restricts the mobility of the system to a submanifold 
A -4 ji(V„+i), implicitly represented by the system 

(3.16) 9 a(t,q\q l ) = 0, a = l,...,n-r, 

with g„ 6 &(j 1 (V n + 1 )) and rank || II = n ~ r - 

The geometric setup developed in Sec. 3.1 is recovered regarding the intrinsic Lagrangian Jzf as the 
pull-back of the extrinsic one. 

In coordinates, the situation is expressed by the equations 

(3.17a) ££{t,q l ,z A ) = C(t,q\if l (t,q % ,z A )), 

(3.17b) g a (t,q\^{t,q\z A )) = 0. 

Let £ denote the product manifold V„ +1 xE n_r , with local coordinates t, q l , X 7 . The natural projection 
£ —4 Vn+i makes £ into a vector bundle over Vn+i and therefore also into a fibre bundle over R. 

We regard £ —> 1R, as the configuration space of a fictitious unconstrained system, merge the functions 
C , g a into a function C := C + X a g a , and adopt the latter as a (singular) Lagrangian over ji(£). 

The content of Eqs. (3.17a,b) is then summarized into the single expression 

(3.18) J?(t,q\z A ) = Cffiq\r(t,q\z A ), A“) V (A 1 , . . . , X n ~ r ) e R” _r . 

17 










From the latter, we derive the relations 


(3.19) 


dT£_ _ dC ch/A dJz? _ d£ dC_ ch/S d& dC dp k 

dt dt dq k dt ’ dq l dq i dq k dq i ’ dz A dq k dz A 


By the arbitrariness of A 1 ,..., A ra r , each equation (3.19) is actually equivalent to n — r + 1 separate 
relations. In particular, the last one splits into 


(3.20) 


d££ dC dp k 

dz A dq k dz A ’ 


dg a dif k 


dq k dz A 


= 0 . 


To any continuous, piecewise differentiable section g : [i 0 ,ti] —» £ let us now associate the action 
integral 


N 


(3.21) 




'jl(v) 


t dt := 


s=l 


C[t,q l 


dq i 

dt 


+ X 7 9<j t, q 


dq i 

dt 


dt 


We have then the following 

Theorem 3.2. The projection £ —->■ V„+i sets up a 1-1 correspondence between extremaloids of the 
functional (3.21) and ordinary extremaloids of the functional (3.1). 


Proof. According to elementary (holonomic) variational calculus, the extremaloids g : q 1 = q l (t ), 
\ a = A a (t) of the functional (3.21) are determined by the Euler-Lagrange equations 


9a(t,q\q z ) = 0 

d dt dt 

dt dq l dq l 


completed by the Erdmann-Weierstrass corner conditions, asserting the continuity of the “momenta” 
and of the “Hamiltonian” H = q l — t along ji(g) ■ 

Setting 7 = 7 Ti ■ /z, po = . Pi = and taking Eqs. (3.17a,b), (3.18), (3.19) into account, 

it is readily seen that the 1 -form p = podt | 7 + Pidq 1 ^ is continuous along 7 and satisfies the conditions 


(3.22a) po + Piip l (t,q\ip l ) = Jf(t,q\%f l ), 


(3.22b) 


dp 0 _ d& dip 1 dpi _ d jgf dif k 
~dt ~ ~dt ~ Pl ~dt ’ ~dt ~ ~d<f ~ Pk ~£kf 


dp k _ d Jzf 

Pk dt A ~dt A 


reproducing the content of Eqs. (3.13). 

This proves that the class p(q) is non-empty, thus ensuring that the section 7 is an ordinary ex- 
tremaloid of the functional (3.1). 

Conversely, if 7 is an ordinary extremaloid, the condition p(y) 7 ^ 0 ensures the existence of at least 
a 1-form p = p 0 dt | 7 +pidq l satisfying Eqs. (3.22a,b). 

From these, taking Eqs. (3.20), the functional independence of the g^s and the Rouche-Capelli the¬ 
orem into account, we conclude that the linear system 


(3.23) 


Pi{t) = 


= +A °(?Sz 


dq l Jhh) 


dga 

dq k 


3 1 ( 7 ) 


= ( d A\ 

dq l ). 


nil) 


admits a unique solution A“ = A“(i), a = 1,... ,n — r. The pair (7 ,p) determines therefore a unique 
section g : [to,t\\ —>• £ , satisfying 7Ti • g = 7. 

On the other hand, through a straightforward pull-back procedure, the 1 form p induces a continuous 
1 -form p = po dt\^ + p t dql along g. 
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On account of Eqs. (3.18), (3.19), (3.22a,b), (3.23), denoting by p the tangent vector field to p, the 
latter satisfies the relations 

(p, p) = po +Pi9| il(/i) = Po +Pi<ptf = = 4h(M)> 

d Po _ /d 2 z?\ _ fdfA _ _ / dC\ fd^\ _ f dt\ 

dt \ dt A P V dt A V dt A \ 9q l / ji ( 7 ) V dt ) - V dt )Hi )’ 

dpi_ _ /dJz?\ _ _ /aSf\ _ fdC_\ (dfA\ _ ( dt\ 

dt V dqi A Pk V dc t i A V dqi A V dqk JMj) V 9q j A V dqi Jhh) ’ 

whence also, making use of the identity {-§fk) u p^ = 0 

dp _ dp 0 dpi i , _ ~ 

jp ~ HA + ~dt dq \j At) + Pi dq \h(rt) ~ dL \HA 

This proves that the class p(p) is non-empty, thus ensuring that the section p is an extremaloid of 
the functional (3.21). □ 

From the proof of Theorem 3.2 one can easily derive the following 

Corollary 3.2. An extremaloid 7 of the functional (3.1) is normal if and only if there exists exactly 
one extremaloid p of the functional (3.21) satisfying 7 = 717 ■ p. 

3.3. Pontryagin’s equations in C(A). The algorithm (3.10a,b,c) involves a set of 2 n + r equations 
for the unknowns q l (t), z A (t ), Pi(t). As such, it has a natural setting in the geometrical environment 
provided by the contact bundle C(A). 

As pointed out in Sec. 2 . 2 , the latter is a vector bundle over A, referred to fibred coordinates 
t , q 1 , z A , pi , and identified with the pull back of the space V*(y n +\) through the commutative diagram 

c(A) y*(Vn+0 

c l A 

A —► V„+l 

The advantage of the environment C(A) comes from the presence of the Liouville 1-form (2.6). By 
means of the latter, the Lagrangian Jzf € A ) may be lifted to a 1-form t>s£ over C(A) according to 
the prescription 7,8 

(3.24) dj? := Jzf dt + 0 = (jzf — pi ip 1 ) dt + pi dq l := — Jt? dt + p t dq l 

The function Jif'it, q l ,Pi, z A ) := Piif % — «if is called the Pontryagin Hamiltonian. The 1 -form (3.24) is 
called the Pontryagin-Poincare-Cartan 1 -form. 

To understand the role of we focus on the fibration C(A) H V ra+ i given by the composite map 
v := 7 r • k. A piecewise differentiable section (X, [i 0 ,i 1 ]) consisting of a finite family of closed arcs 

X (s) : [a s _i,a s ] ->• C(A ), s=l,...,N, t 0 = a 0 < ai < • • ■ < a N = ti 

is called v-continuous if and only if it projects onto a continuous, piecewise differentiable section 
v ■ X : [to, ti] —t V„+i • 

A deformation X{ = {(X^, [a s _i(£), a s (£)]) } is called u-continuous if and only if all sections Xj are 
u-continuous. The necessary and sufficient condition for this to happen is the validity of the matching 
conditions 

(3.25) vX ( f(a s (^)) = vX[ s + 1 \a s (0) V |£| < £, 0 = 1 ,..., N - 1 . 

A deformation X% is said to preserve the endpoints of v • X if and only if the projection v • Xg does 
so. A vector field along X tangent to the orbits of a -^-continuous deformation is called an infinitesimal 
deformation. 

^ A deeper insight into the geometrical meaning of the 1-form (3.24) comes from the study of the gauge-theoretical 
structure of the calculus of variations, as developed in [39]. 

^ As usual, we use the same symbol for covariant objects in A and for their pull-back in C(A). 
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It is worth noticing that the previous definitions do not require the admissibility of the sections v • X 
or v ■ In this respect, the only condition needed in order for a vector field X 1 + r‘ 4 (^r) x + 

n*(g§T) x t° represent an infinitesimal deformation of X is the consistency with the matching conditions 
(3.25), summarized into the jump relations 


(3.26) 


X l + a, 


dq i 

dt 


= 0 


s = 1,..., iV — 1 


with a s = (^fr)j _ 0 and with no a-priori relationship between and , 0 ( x - 

By means of we define an action integral X[X\ over C(A), assigning to each u-continuous section 
X : q l = q l (t ), = z A (t ), pi = pi(t) the real number 


(3.27) 


X{X\ := / $se = 


/to 


p^~^) dt = Y. 


N 


— 1 J a s - 1 \ 




|x( s 


dt. 


For any u-continuous deformations X^ preserving the endpoints of v -X, denoting by + 

+ IIi(^r) the corresponding infinitesimal deformation, we have then the relation 


dZ[Xff\ 


f 

( dq i 


f d Pi 


\x< 6 *t 4 

df 

£=o 

Jt 0 

[\ dt 

fa) 

\ dt 

dq 1 j 

dz A 


dt 


N-1 

£ 


Pi ( X 2 + a s ) - a s Jf? 


dX I 


From this, taking Eqs. (3.26) into account, we conclude that the vanishing of 
is mathematically equivalent to the system 

dq i dJ? d Pi dJ4? dJA 

dt 


for arbitrary Xt 


(3.28a) 


dt 


dq i 


dz A 


= 0, 


s = 1 ,... ,7V — 1 . 


dpi 

completed with the continuity conditions 

(3.28b) [ Pi \ as = [je] at =0 

In view of the definition of the Pontryagin Hamiltonian S4?, Eqs. (3.28 a,b) are easily seen to reproduce 
the content of Eqs. (3.10a,b,c), (3.11) 9 . 

As far as the ordinary extremaloids are concerned, the original (constrained) variational problem in 
the configuration space is therefore equivalent to a free variational problem in the contact manifold, in 
full agreement with Pontryagin’s maximum principle. 

As a further comment on Eqs. (3.28), let us digress on the special situation determined by the ansatz 
££ — 0. In the stated circumstance, the functional 


(3.29) 


Z 0 [X] := / 0 = 


I to 




is an intrinsic attribute of the manifold C(A), entirely determined by the Liouvillc 1 form (2.6). Its role 
is clarified by the following 

Proposition 3.1. Let 7 : [to,ti] —>• V n _|_i denote any continuous, piecewise differentiable section. Then: 

a) 7 is admissible if and only if the functional (3.29) admits at least one extremaloid X projecting 
onto 7 , i.e. satisfying v ■ X = 7 ; 

b) for any admissible 7 , the extremaloids of Iq projecting onto 7 are in 1-1 correspondence with 
the elements of the annihilator (Y( 2 H))°. 


^Actually, exactly as it happened in Sec. 3.1, the condition [ q' j =0 does not arise from the stationarity requirement, 
but is implicit in the definition of X. 
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Proof. For ££ = 0, 

Eqs. (3.28a,b) reduce to 

(3.30a) 

aJ 

II 

‘A. 

(3.30b) 

dpi dif k 

~dt + ~dT Vk ~ 

(3.30c) 

dib 1 

(3.30d) 

[Pi\a e = [PiP] 


Eq. (3.30a) is the admissibility requirement for the section v ■ X. Accordingly, if an extremaloid X of 
the functional (3.29) projects onto 7 , its image ( ■ X under the map ( : C(A) —» A coincides with the 
lift 7 : [t 0 , £ 1 ] — A. 

For any admissible 7 , the extremaloids X projecting onto 7 are therefore in 1-1 correspondence with 
the solutions pi(t) of the homogeneous system (3.30b,c,d), with the functions q l (t), z A (t ) regarded as 
given. 

Assertions a) and b) are then straightforward consequences of the fact that, according to Proposition 
2.4, Eqs. (3.30b,c,d) are precisely the conditions required in order for a virtual 1 forms Pi(t) 6q z | 7 to 
belong to the annihilator (T ( 2 U)) 0 . □ 


According to Proposition 3.1, a section 7 : [£ 0 ,£i] —> V n +i describes a normal evolution of the system 
if and only if the functional (3.29) admits exactly one extremaloid projecting onto 7 , namely the one 
corresponding to the trivial solution Pi(t) = 0. If the extremaloids projecting onto 7 are more than one, 
7 represents an abnormal evolution; if no such extremaloid exists, 7 is not admissible. 

Returning to the action integral (3.27) we can now state 


Corollary 3.3. The totality of extremaloids of the functional (3.27) projecting onto a section 
7 : [£o,£i] —>• V ra +i is an affine space, modelled on the vector space formed by the extremaloids of the 
functional (3.29) projecting onto 7 . 


The proof, entirely straightforward, is left to the reader. 

The previous arguments provide a restatement of Theorem 3.1 in the environment C(_4). As already 
pointed out, the projection algorithm X —» v-X does not reproduce all the extremaloids of the functional 
(3.1), but only the ordinary ones. 

The missing ones may be obtained determining the abnormal evolutions by means of Proposition 3.1, 
finding out which ones have an exceptional character, and directly testing the validity of Eq. (3.2a) on 
each of them. 


3.4. Hamiltonian formulation. As pointed out in Sec. 3.3, all ordinary extremaloids of the functional 
(3.1) are projections of extremaloids of the functional (3.27). Let us discuss the implications of this fact. 

To this end, temporarily leaving aside all aspects related to the presence of corners, we observe that a 
curve X in C{A) is an extremal for the functional (3.27) if and only if its tangent vector field Z := X* ( 
satisfies the properties 


(3.31) 


{Z,dt) = 1, ZA d'dse = 0. 


On account of Eq. (3.24), at each q £ C(A), a necessary and sufficient condition for the existence of 
at least one vector Z £ T^{C{A)) satisfying Eqs. (3.31) is the validity of the relations 


(3.32a) (|5) = °’ A = 1 ,... ,r. 

Points at which Eqs. (3.31) admit a unique solution Z will be called regular points for the functional 
(3.27). In coordinates, the regularity requirement is expressed by the condition 


det 


/ d \ 
\dz A dz B J ? 


^0. 


(3.32b) 
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In view of Eq. (3.32b), in a neighborhood of each regular point Eqs. (3.32a) may be solved for the 
z A ’s, giving rise to a representation of the form 

(3.33) = z A {t,q 1 ,...,q n ,p 1 ,...,p n ). 


The regular points form therefore a (2 n+ 1 (-dimensional submanifold S A C{A), locally diffeomorphic 
to the manifold V*(V n +i)- 

Inserting Eqs. (3.33) into the first pair of relations (3.28a) gives rise to a system of ordinary differential 
equations in normal form for the unknowns q z (t),pi(t) . The algorithm is readily implemented denoting 
by TL := the pull-back of the Pontryagin Hamiltonian, expressed in coordinates as 


<f, z A {t, q\ Pi), Pi) = Pk ip k (t, q\ z A (t, q\pi))-^{t, q\ z A (t, q\pf)) 


In view of Eqs. (3.32a) we have the identifications 

dU _ dje_ dJfdz^ _ i 

dpi dpt + J$z A dpi 

dU ddtf dJfdz A dip k dS£ 

dq i ~ dq i + /6z A dq i ~ Pk dq i dq * ' 

On account of these, the first pair of equations (3.28a) takes the form 


(3.34a) 

(3.34b) 


dq 1 dU 

dt dpi ’ 

dpt _ d'H 
~dt -~~dq 1 ' 


The original constrained lagrangian variational problem is thus converted into a free hamiltonian 
problem on the submanifold j : S —> C(A) ■ 

A u-continuous extremaloid of the functional (3.27) consisting of a finite family of closed arcs : 
[a s _i,a s ] —>• C(A), each contained in (a connected component of) the submanifold S will be called a 
regular extremaloid. 

Singular extremaloids, partly, or even totally lying outside S may also exist. In fact, while Eq. (3.32 
a) is part of the system (3.28a,b), and must therefore be satisfied by any extremaloid, the requirement 
(3.32b) has only to do with the well-posedness of the Cauchy problem for the subsystem (3.28a). 

On the other hand, by construction, the Hamilton equations (3.34a,b) determines only the regular 
extremaloids. The singular ones, when present, have therefore to be dealt with directly, looking for 
solutions of Eqs. (3.28a,b) not arising from a well-posed Cauchy problem. In principle, this could be 
done extending to the non- holonomic context the concepts and methods commonly adopted in the study 
of singular Lagrangians [40]. The argument is beyond the purposes of the present work, and will not be 
pursued. 


Appendix A. Finite deformations with fixed end points: an existence theorem 

(i) Given an admissible, piecewise differentiable section 7 : [to,H] —>• V n +i, a crucial question 
is establishing under what circumstances every admissible infinitesimal deformation vanishing at the 
endpoints of 7 is tangent to an admissible finite deformation 7 ^ with fixed endpoints. The following 
preliminaries help simplifying the discussion. 

Proposition A.l. Let 7 : (c, d) —> A be the lift of an admissible differentiable section 7 : (c, d) —> V„+i . 
Then, for any closed interval [a, 6] C (c, d) there exists a fibred local chart (U , k), k = (t, q 1 ,..., q n , z 1 ,..., z r ) 
satisfying the properties 

(A.la) • 7 (t)GU V f e [a, b\; 

(A.lb) • the intersection 7 ((c,d)) 0 U coincides with the curve q l = z A = 0; 

(A.lc) . ^(7 (t)) = (§^)- (t) = 0 V 7 W G U. 
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Proof. The existence of fibred local charts (U, k), k = (t, q 1 , ■ ■ ■ q n , z 1 , ■ ■ ■ z r ) satisfying Eqs. (A.la, b) 
and the first relation (A.lc) is entirely straightforward. 

Choose any such chart, and denote by q 1 = tf'ft, q l , z A ) the corresponding representation of the 
imbedding A — > j i(V n +i) • 

Under an arbitrary transformation q l = qf z A = z A we have then the transformation laws 


i > 1 = q° + a 1 ^ 3 , 


dfp 

dq k 


da^ 

dt 


i 


dqi 


(OV 


Therefore, if the matrix a l .j (t) obeys the transport law 


da l : 

dt 


+ a 1 


w)~ 


o, 


7 b) 


the coordinates t,q‘,z A satisfy all stated requirements. 


□ 


Every local chart (U,k) consistent with Eqs. (A.la,b,c) is said to be adapted to the arc ( 7 , [a, 6 ]). 
For later use we observe that this property is stable under arbitrary transformations of the form 

(A.2) <f = <f(g\ • • • ,q n ), z A = z A (t,q 1 ,...,q n ,z 1 ,...,z r ), 

with ( 0 ,..., 0 ) = z A (t, 0 ,..., 0 ,..., 0 ) = 0 . 

Corollary A. 1. Let 7 = {( 7 ^, [a s _i,a s ]), s = l,...,7V} be the lift of an admissible piecewise differen¬ 
tiable section ( 7 , [to,ii]). Then, there exist fibred local charts (U s ,k s ), k s = (t, q A ,..., gA , z A, ..., z/\) 

adapted to the arcs 7 ^) such that, in each intersection tt(U s ) f~l 7 r(/ 7 s+ i), the coordinate transformations 
q{ s + 1 ) = 9 (g + 1 ) {t, 9 (g), • • • j 9 (")) satisfy the conditions 


(A.3) 


dq ( \ 


(» + l) 


d( l(s) ' 7(a s ) 


= *5 


dqt 


( 8 + 1 ) 

dt 


7 K) 


V’XtO 


Proof. The conclusion follows at once from Proposition A.l, observing that the freedom in the choice of 
the the adapted coordinates summarized into Eqs. (A.2) leaves full control on the values of the Jacobians 

( ( )' P art i c ular, i 11 the intersection n(U s -i Cl U s ), the arcs 7 ^, 7 ( s + 1 ) are respectively 

described by the equations 

qls + i)(7 (s) ) = 9(s + i)(t,g( 1 s) (7 (s) ),---,9(”)(7 (s) )) = «(» + i)(t,°,...,0), 

9(U i)(7 (s + 1) ) = 0. 

In the coordinate system t, q (s + i)’ the jump of the tangent vector at the corner 7 (a s ) is therefore 
given by 

'^( S + 1)(7 (S + 1) ) dq( s + i)( 7 (s) )\ _ (dql + i) 


^(7) 


dt 


dt 


7(“«) 


dt 


l( a s) 


D 


Every family of local charts { (U s , k s ) , s = 1 ,..., A} satisfying the requirements of Corollary A.l is 
said to be adapted to the section 7 . 

(ii) Assigning an adapted family of local charts singles out a distinguished class of controls 
< 7 ^ : n(U s ) —> U s , described in coordinates as 

(A.4) ff (s) : zfa(t, qfa ,..., qfc) = 0. 

By construction, each such <7 M satisfies a ^ -7 = 7 , i.e. it contains the section 7 in the sense 

described in Sec. 2.1. 

For each s = 1, ...,1V, the restriction of the tangent map eri^ to the vertical bundle along 7 7) 
determines an infinitesimal control h^ : U( 7 ^ s -*) —> expressed in coordinates as 
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hi A {t) = 0 . 




d lls) 


d s ) (t) 


= cri^ 




A s ht) 


d < S ) 


d s )(t) 


In view of Eqs. (2.24b), (2.25a,b), (A.lc), the absolute time derivative associated with satisfies 


(A.5) 


D 


D 


= 0, 


s = l,...,N. 


^V^,)A(.) (i ) Dt I'M*) 

Noting that, according to Eq. (A.3), the fields — and therefore also the virtual 1-forms 

Sq l j 7 (s) — match continuously at the corners, we conclude that the sections e(j) : [to,ti] —> V( 7 ), eW : 
[to,ti] —► V*("f) respectively defined by 

d 


(A. 6 ) 


e(i)(0 = 


e(*)(t) = (Sg* 7W(t) , 9 = 1 ,..., 


form a dual bases for the spaces 14 , 14 * of h-transported helds along 7 . 


AT 


(iii) Let us now come to the main question. Let 7 := {( 7 ^, [a s _i,a s ]) } be an admissible piecewise 
differentiable section, {(f7 s ,fc s )} a family of local charts adapted to the lift 7 , and {e(j)}, {e^} the 
corresponding A-transported bases. 

We recall that, with the notation of Sec. 2.7, the most general infinitesimal deformation X of 7 is 
determined, up to initial data, by an element ([/, a) € 2U, namely by a vertical vector held U along 7 and 
by a collection of real numbers a = (a ±,..., oat-i) . In particular, a necessary and sufficient condition 
for X to vanish at the endpoints of 7 is expressed by the requirement T (U,a) = 0 which, in adapted 
coordinates, reads 


(A.7) 




N—l 


dt - ^a s [y ( 7 )] 


= 0 . 


To inquire whether a given infinitesimal deformation vanishing at the endpoints of 7 is tangent to a 
finite deformation with fixed endpoints we introduce an auxiliary tool, namely a positive-definite scalar 
product in V (A). 

From a tensorial viewpoint, this means assigning a contravariant held over A , described in hbred 
coordinates as G = G AB gfa-® -jpw, with the components G AB identifying a symmetric, positive-dehnite 
matrix, inverse of the matrix Gab formed by the scalar products (gdr, ~ipw) ■ 

Given any {U,a) G 2H, we shall now construct a family of hnite deformations 7 (g,i/) of the original 
section 7 , tangent to the inhnitesimal deformation determined by (G,a) and depending on n auxiliary 
parameters, identihed with the components of a /i-transported 1-form v £ V£ . To this end, we introduce 

• a family of functions 


(A.8a) a s (£,v) := a s + a s £ - \ a s 2 £ 2 [^ l {l)\ ae v% , s = 1,..., TV - 1 , 

expressing the deformation of the partition of the interval [to,ti] into subintervals [a s _i,a s ]. For 
notational convenience, the quantities (A.8a) are completed by the constant functions 
a o(£j v ) = t 0 , a N {^u) = 4; 

• a family of deformations : n(U s ) —> U s of the controls (A.4), described in adapted coordi¬ 
nates as 10 


(A. 8 b) 


4 =^f.)(*) + (s) := zfa{£,v,t). 


A straightforward argument shows that for any bounded open subset A C V£ there exists m > 0 
such that each image (tv(U s )) is contained in U s for all v G A, |£| < m. 


■^Eq. (A.8b) is strictly coordinate-dependent. As such, it has no invariant geometrical meaning, but is merely a tool 
for the subsequent construction of a family of finite deformations. 
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Theorem A.l. Let 7 be an admissible piecewise differentiable evolution. Given any {U,a) £ 2T, define 
the functions as^jis) and the sections cr^ a K as above. Then, for any open bounded subset A C Vf 
there exist an e > 0 and a family 7(£, v ) = { ( 7 (|y) ’ [ a s-i(£; a s(£> z/ )]) } of piecewise differentiable 

admissible sections defined for |£| < e, v £ A and satisfying the properties 

a) 7(0,i/) (t) = l(t) Vv; 

b ) 7 («,i/)(*o) = l(t 0 ) V£, zz; 

c ) = T {a s {£,,v)) V£, v , Vs = 1,..., N - 1 

d) each arc 7u\-) is contained in the control cr^ s ^, namely it satisfies the condition = 

(a) (s) 

‘’(UrUV 


Proof. As pointed out in Sec. 2.1, each control (A. 8 b) determines a velocity field over n(U s ), depending 
parametrically on £ and 1 /. The latter may be viewed as a vector field in the 

product manifold (— m,m) x A x 7 t(U s ), with components ZX(£,v,t,qX) = if l (t, qA,zA(^,ts,t)). 
,( s ) 

Let us denote by Q ; the local 1 parameter group of diffeomorphisms generated by Z/ s ), and by x s 
the corner 7 (a s ). 

Then, on account of Eq. (A.lc), for any v* £ A the orbit of through the point (0, v*,x s -\) 
coincides with the coordinate line q 1 = 0 , £ = 0 , v = v*. and is therefore defined for all t in an open 
interval (b s _i,b s ) D [a s _i,a s ]. 

Taking the compactness of A into account, a standard result in classical Analysis [32, 41] ensures the 
existence of a positive e and of open neighborhoods W s -± 9 27 - 1 , s = 1,..., A, such that each map 
is well defined on (—e, e) x A x W s -± for all t in the closed interval [a s _i(£, v), a s (£, v)\ C (b s _i, b s ). 

From this, denoting by E s the slice t = a s (£, z/) in (—m, m) x A x V n +i, we conclude that the group 
(s) 

Q maps the intersection W s -1 fl E s _i into E s . Without loss of generality we may always arrange that 
the image of each W s _i fl E s _i is contained in W s fl E s , s = 1,..., N. 

The rest is straightforward: for each |£| < e, v £ A, consider the sequence of closed arcs 7^1) : 
[a s _i(£, v), a s (£, v)] —> tt(U s ) defined inductively by 

7= Ct ( 1 ) (^ z/ >7(*o)) t £[t 0 ,ai{€,u)], 

7(^{ 1} W = Ct (S + 1) (^^,7 <£l)( a s(t,v))) t£[a s (^),d s+ 1(0] ■ 

The collection 7 ^ ^ := {( 7 ^^, [a s _i(£, z/), a s (£, z/)]), s = 1 ,...,1V} is then easily recognized to 
fulfil all stated requirements. □ 


(s') 

In adapted coordinates, each arc 7 ^ is represented in the form 

(A.9) q{ s) = , a s _i (£,v) <t<a s ((;,v), 

with the functions 7 >f s \ satisfying the equations 

is) 


(A.10a) 


dt 


r(t, (ufa t ) + w(g ab |£) i(<) J 


as well as the matching conditions 

(A.10b) <p{ s + 1) (f„v,a s {t,v)) = gp, + i)(a s (^zz),^ ( l s) (^i/,a s (^j/))), 

< 7 (s + 1 ) = 9(8 + i)(^i 9(s) i • ■ • i 9(s)) denoting the transformation between adapted coordinates in the in¬ 
tersection n(U s n U s +i). 

From this, taking Eqs. (A.lc) into account, it is easily seen that, for any is* £ A, the 1-parameter 
family of sections 7 ({,j/*) is a deformation of 7 , tangent to the infinitesimal deformation X determined 
by the vector ( U , a) £ W. 

/ dip / \ ,£) 

We have therefore the identifications ( —- 

(2.33). 


= X (s)(i) i completed with the jump relations 
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After these preliminaries, let us now focus on two facts: 

• on account of Theorem A.l, given any bounded open subset A C V f *, the correspondence 
(£, u) —>• 7(£,i,)(ti) determines a differentiable map X of the cartesian product (— e,e) x A into 
the slice t = t\ in V n +i, with values in a neighborhood of the point 7(^1); 

• given any differentiable curve v = v(£) in A, the 1-parameter family of sections 7(f,i/(f ))(£)> 

|£| < e, t £ is a deformation of 7, leaving the first endpoint 7(^0) fixed. 

Both assertions are entirely obvious; in adapted coordinates, the map X is represented in the form 
(A. 11) <l(N)( x (^v)) = ^(ao(£>Mi) := 

Exactly as above, it may be seen that, independently of the choice of the function ^(£), the deformation 
7(£>"(£)) always tangent to the infinitesimal deformation X determined by the vector ( U,a ) £ 2U. 

From this, taking the relations X*(0,^) = 0, (=^-) tQ = X l (ti) into account and recalling Taylor’s 
theorem, we conclude that, whenever the condition X(t\) = 0 holds true, i.e. whenever the vector ( U,a ) 
belongs to ker(T), the left-hand side of Eq. (A.11) may be cast into the form 

(A.12) X j = i Z 2 0%v), 

with regular at £ = 0. 

In this way, the original problem is reduced to establishing under what circumstances the valid¬ 
ity of Eqs. (A.7) is sufficient to ensure the existence of an s' > 0 and of a curve i/(£) satisfying 
*(£, KO) =7(*i) V|£| < s'. 

In adapted coordinates, on account of Eq. (A. 12), the answer relies on the solvability of the equations 
(A.13) 0 % (£,vi,...,v n ) = 0 i = 1,... ,n 


for the unknowns 17 in a neighborhood of £ = 0. 

To examine this aspect we notice that, in each adapted chart, Eqs. (A.10a) imply the transport law 


d_ 

dt 


d 2 V( s ) 

d£ 2 


€=o 


/ d 2 ^ 

\dq k dq r J -(s) 


X k X r 


d 2 ip l 
dq k dz A J -{ s 


X k U A + 


/ d 2 ^ 1 
V dz A dz B 



U A U B 



/ d^ 1 \ 


M G 


-<AB 


dip k \ 
dz B )^(s) ’ 


the cancelation arising from Eq. (A.lc). 

In a similar way, from the matching conditions (A. 10b), recalling Eqs. (A.3), (A.8a) and evaluating 
everything at £ = 0, we get the relations 


d 2 'A 


^(« + P 

<9£ 2 


€=o 


d 2 vl) 

d£ 2 


C=o- 


r ) 2 (-j i 

2 ° q (s+l) 


r) 2 n i 

+ !) y h X k o n . 

— X (s) X (s) - 


dt 2 

dX 


2a, 




q ( s + 1) \r k 


dt dq, k X ^ + 


'W 


dq { h s) dq k g) 


(s + 1) dX (s) 


expressing the jumps 





LI a? 2 4=o 

\ dt 2 , 

S=oJ 


dt dt 

at the corners. 




^(7 ) ^(t) 


Collecting all results and recalling Eqs. (A.11), (A.12) we obtain the expression 
'd 2 X^ 

h =0 


l«=o 


dt 2 


N -1 




(A.14) = b l 

with b l £ R independent of v. We can therefore state 


^*(7) ^(7) 


\Vk 


Vk 
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Proposition A.2. Let 7 : [fo,ti] —> Vn+i be a continuous, piecewise differentiable, admissible section. 
Then, if the matrix 


(A.15) 



ab df_ d£\ 
dz A dz B 


dt 


N-l 

E«. 2 

S— 1 


^(7) ^*( 7 ) 


is non-singular, every infinitesimal deformation of 7 vanishing at the endpoints is tangent to a finite 
deformation with fixed endpoints. 


Proof. The conclusion follows at once from the fact that, on account of Eq. (A. 14), the non-singularity 
of the matrix (A.15) ensures the solvability of Eq. (A.13) in a neighborhood of £ = 0 , Mu. □ 


Proposition A.2 may be rephrased in the language of Sec. 2.7 observing that, whenever the section 7 
is abnormal, the matrix (A.15) is necessarily singular. 

In the stated circumstance, in fact, Proposition 2.4 and Eq. (A. 6 ) imply the existence of at least one 
non-zero virtual 1 -form pi 5q l | 7 with constant components pi obeying the relations 

(A.i6) «(§?)„„ = 0 ' = » 

and therefore automatically satisfying pi S = 0 . 

More specifically, denoting by p the abnormality index of 7 , we have the following 


Theorem A.2. The matrix (A.15) has rank n — p. 

Proof. By definition, the index p coincides with the dimension of the annihilator (T(2H))° C V)* which, 
in turn, is identical to the dimension of the space of constant solutions of Eqs. (A. 16). 

On the other hand, by Eq. (A.15), the matrix S ® J ' is positive semidefmite. Its kernel is therefore identi¬ 
cal to the totality of zeroes of the quadratic form piPj , i.e. to the totality of n-tuples (p \,..., p n ) £ R” 
satisfying the relation 


0 = 


'to 


G 


AB 


difi 1 
dz A dz B J 


N -1 


d^ k \ J ^ 

707 dt + J 2 


U a 


ril) 1 ^( 7 ) 


Pi Pk — 


'to 


AB or d^ k \ 

G Pi O .4 Pk " ' 


N-l 


dz A ^ dz B ) fi 


dt + 


a s I Pi 


s=l 


rift) 


Due to the positive definiteness of G AB {t ), the last condition is equivalent to Eqs. (A. 16). This proves 


dim ( ker(S'* :? )) = p =>- rank (S*- 7 ') = n — p. 


□ 


Proposition A.2 and Theorem A.2 show that the normal evolutions form a subset of the ordinary ones, 
thus establishing Proposition 2.5 of Sec. 2.7. 

Along the same lines, a useful generalization is provided by the following 


Theorem A.3. Let p (> 0) denote the abnormality index of 7 . Then a sufficient condition for every 
infinitesimal deformation vanishing at the endpoints of 7 to be tangent to a finite deformation with fixed 
endpoints is the existence of an (n—p) -dimensional submanifold S C V„+i such that every deformation 
75 leaving 7 (^ 0 ) fixed satisfies 7 {(ti) £ S for £ sufficiently small. 

Proof. Using the freedom expressed by Eq. (A.2), we choose the adapted coordinates in such a way that, 
in a neighborhood of 7 (U), the submanifold S has local equation q j^ 1 = ■ • ■ = q = 0 n . 

Given an infinitesimal deformation X generated by an element ( U,a ) £ 2H, we then proceed ex¬ 
actly as above, ending up with a finite deformation 7 ({ jt /) tangent to X, described in coordinates as 
qfs) ~ T( s ) (£; v i t) > a s-i(£,z / ) <t< a s (£,v), with the functions ip^i^,v,t) satisfying all conditions of 
Theorem A.l. 


Namely, we choose the last chart (UN,k]y) consistently with the stated condition, ant then proceed backwards, 
requiring the validity of Eqs. (A.4) at the corners. 
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Once again, we focus on the “end-point map” X(£,p) = 7(g,i/)(ii) and on the fact that, under the 
assumption X{t\) = 0, the functions X*(£, p) = ip, l N X£,u,ti) factorize into the product 

* = \eo\^), 

with 6 l | £_ 0 expressed as a linear polynomial 


6 l 


€=o 


d 2 r 


= b l + S i:, v j . 


e=o 


Actually, in the present case, the condition X a (£, v) = 0, a = p + 1,... ,n implies 0“(£, v) = 0 
whence, in particular, b a = S aj = 0. 

At the same time, Theorem A.2 ensures the equality between ranks '*- 7 and the abnormality index p 
of 7 . 

In the case in study, the matrix S^ is therefore necessarily of the form 

/ S AB 0^ 


S 13 = 


0 0 


A,B = 1,... ,p, 


with det S AB ± 0. 

The rest is now straightforward: in order to establish the existence of a finite deformation with fixed 

endpoints tangent to X we have to verify that the equations 0*(£, v) = 0 admit at least one solution 

v = p(£) in a neighborhood of £ = 0. And indeed, no matter how we choose the functions i/ a (£), 

a = p+ 1 ..., n, the n — p equations 0 “(£, v) = 0 are identically satisfied, while the remaining ones form 

a system of p equations for the unknowns i/±,... ,v p , whose solvability is ensured by the non singularity 

of the Jacobian |p— I = S AB . □ 

dv B I {=o 


Appendix B. Comments on normality 

The following arguments help clarifying some aspects of the concept of normality discussed in 
Sec. 2 . 7 . 

Let 7 = {(7W, [a s _i,a s ])} be a piecewise differentiable evolution. According to Proposition 2 . 4 , if 
at least one arc 7^ is normal, 7 is necessarily normal. 

More generally, an evolution may happen to be normal even when all its arcs 7W are abnormal. 
Examples in this sense are: 

• V n +i = R x E2, referred to coordinates t,x,y. Constraint: x 2 + y 2 = v 2 . Imbedding A —> ji(V n +i) 
expressed in coordinates as x = v cos z, y = v sin z■ Piecewise differentiable evolution 7 consisting of 
two arcs: 

7^: x = 0 , y = vt t o <t <0 

7 ^ : x = vt, y = 0 0 < t < t\ 

Eq. ( 2 . 42 a) admits /i-transported solutions p^ = aSy | 7 , p^ = fiSx^ (a, (3 G R) respectively along 7^) 
and 7^: both arcs are therefore abnormal. Nevertheless 7 is normal, since no pair p( 1 \p( 2 '> matches 
into a continuous non-zero virtual 1 -form along 7. 

• V„+i = R x E 2 , referred to coordinates t,x,y. Constraint: v 3 x = ( y 2 — a 2 t 2 ) 2 . Imbedding 
A —> ji(V„+i) expressed in coordinates as x = v ~ 3 [z 2 — a 2 t 2 ) 2 , y = z. Piecewise differentiable evolution 
7 consisting of two arcs: 


(!) . 

x = 0 , 

y 

= \a(t 2 -t* 2 ) 

t 0 <t <f 


„ 4 




( 2 ) . 

x = — (t 5 -C 5 ), 
bv 6 

y 

= 0 

t* <t<t 


(t* ^ 0). Eq. (2.42a) admits /i-transported solutions of the form p = aSx i 7 along the whole of 7 . Both 
arcs 7 ^-*, " 1 ^ are therefore abnormal. Nevertheless, 7 is normal, since no solution satisfies condition 
(2.42b). 
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• By definition, local normality implies normality. The converse is generally false, as shown by the 
previous examples. A further example, not involving the presence of corners, is the following: let the 
imbedding A ——>■ j i(V„+i) be locally described by the equations 

( q A = z A A = 1 ,... ,n — 1 

l q n =f(t)z 1 

with f(t) = exp(—i/t 2 ) for t < 0 and f(t) = 0 for t ^ 0. 

Along any admissible section 7 : [to,ti] V n + 1 , the condition of ^.-transport and Eqs. (2.42a) are 
summarized into the relations 

(B.l) —— = 0, Ai + A„/(f) = 0 , A 2 = • • • = A„_i = 0 . 

at 

In particular, if to < 0 < t\ we conclude that: 

o 7 is normal, since Eqs. (B.l) do not admit any non-zero solution in [to,ti]; 
o 7 is not locally normal, since Eqs. (B.l) admit the solution Xa = 0, A ra = const, in any subinterval 
[M] C [0, . 
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